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Abstract: Various aspects of time-dependent processes are studied within the large N 
approximation of 0(N) vector models in three dimensions. These include the rolling of 
fields, the tunneling and decay of vacua. We present an exact solution for the quantum 
conformal case and find a solution for more general potentials when the total change of 
the value of the field is small. Characteristic times are found to be shorter when the time 
dependence of the field is taken into account in constructing the exact large N effective 
potentials. We show that the different approximations yield the same answers in the regions 
of the overlap of the validity. A numerical solution of this potential reveals a tunneling in 
which the bubble that separates the true vacuum from the false one is thick. 
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1. Introduction 

Evolutionary processes were studied in particle physics during the whole time of its own 
evolution. However, in recent years the study of time-dependent processes has played 
an increasing role in cosmology as well as Quantum Field Theory. The spot light has 
shifted away from the study of vacua and perturbations around it. Metastable vacua and 
their properties are investigated in detail. In many cases an understanding of the time- 
dependent processes is important. In addition to the interesting theoretical issues involved, 
the problem emerged in cosmology in general, and in the cosmology which deals with a 
notion of a landscape of vacua in particular. The time-dependent processes determine both 
the tunneling to a different potential well and the rolling of a field down to the region of a 
minimum of a potential. Similar physical processes appear also in some field theory models 
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of possible phenomenological interest, like for instance in the supersymmetric models with 
long-living metastable vacua in which SUSY is broken [P . Also in the Standard Model 
and beyond it such a tunneling followed by a rolling towards the minimum can occur. 

In most cases a detailed study of these processes is difficult. Systems accounting for 
gravity and ignoring it were studied. Various approximations have been used, such as the 
semiclassical and the thin wall approximations ||2|, ^ . 

Understanding the exact features of such time-dependent processes is important for 
theory and may eventually have practical uses. 

For this purpose we conduct a study in the case of large vector models in three 
space-time dimensions. 

Traditionally, quantum field theories with N dynamical variables, — > oo, have served 
mostly as study grounds for extending one's intuition in handling basic problems of quan- 
tum field theories |5|. Approximate results obtained at large N possess many of the proper- 
ties believed to be true in an exact solution. Such processes have been studied in the large 
A^ approximation of the vector model in d = 4. Here we obtain new results and highlight 
additional aspects for d = 3, which includes a quantum conformal case. 

In this paper we investigate various time-dependent processes in this model. In section 
|2| we review the rich phase structure of the three dimensional vector models and the methods 
of solving them. We next calculate in section ^ an exact large N effective action which 
allows treating time-dependent solutions. In the derivation we take into account time 
variations also in the potential terms. We derive the general equations which govern an 
evolution of the system and solve them in two special cases. In the first case (section 
we consider a rolling of the field from the top of a potential (j)^. The corresponding 
theory is also quantum conformal [^, and it makes a solution accessible. We show that the 
field continues to escape to infinity in a finite time. A time is shorter, however, than that 
obtained in an approximation where the time dependent effects in the potential are not 
accounted for. This situation is of interest in studying possible holographic descriptions of 
Big Crunch singularities Q . 

In the second case (section ^ we consider a rolling in a general potential in an approxi- 
mation of small changes in fields along the evolution. We derive a characteristic time of run 
away from a potential maximum and a frequency of small oscillations around a minimum. 
The frequency of oscillations is bigger than the corresponding frequency obtained for an 
effective potential derived assuming static solutions only. 

We also study a possibility of a tunneling in the system (section ^). We show that a 
solution in this case looks like an expanding bubble of a complicated shape different than 
that resulting in the thin-wall approximation |Q] in the semiclassical case. We illustrate 
this fact by a numerical example. 

We end with a set of appendices. 

2. Scalar model in the large N limit - A brief review 

Let us consider an 0{N) symmetric Euclidean action for an A - component scalar field (j) 
in three space-time dimensions 
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s 



(2.1) 



The potential has a Taylor expansion of the form 



n=l 



(2.2) 



with g2n kept fixed as ^ co (in the large N hmit any such potential is renormalizable 
B). We describe a way to find the generating functional of this model. It is given by 



Z 



J 



Dcpexp 



5(0) - / J{x) ■ (l3{x)<fx 



Inserting 

1 ~ j Dp5{^ - iVp) ~ J DpDX 
and integrating over (j) one obtains 



e J 2 



J 



DpDX e 



-NS,ff{p,X) y J{x){-n+iX)-^Jiy)d^xd3y 



where 



Seffip, ^) = \j \'^U{p) - iXp] d?x + ^Tr In (-□ + i\) . 



(2.3) 



(2.4) 



(2.5) 



(2.6) 



When is large, the last form suggests using the saddle point method to calculate the 
integral. The two saddle point equations, obtained by varying the auxiliary fields p and A 



are 



2U'{p) =iX, p = tr- 



1 



/iX 



m 
47r 



— (2 7) 

-n + iX in 47r ' ^ ' ^ 

where m will assume a role of a mass. Moreover, we have used the dimensional regular- 

ization procedure in order to define the divergent loop. One can also define an effective 

potential whose minimum fixes a value of the mass 



U, 



eff 



N 



n=l 



(2.8) 



In the above analysis and in what follows we take into account the fact that there is no 
spontaneous breaking of the 0{N) symmetry. The complete vacuum energy ^ which 
accounts for 0{N) broken phase is given by 



eff 



N 



247r 



n=l 



2n 



$1 
N 



m 



(2.9) 



""^We use here the definition Tr = J dPxtr 
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where is the would be expectation value of (j). 

When one considers only relevant and marginal terms up to (j)^, then the potential is 
stabilized for < ge < (47r)^. It is unstable for qq < 0, the instability implied for qq > Qc 
is discussed in the next section. If only is present the theory is conformal and possesses 
two 0{N) invariant phases, one with g^ < (47r)^ and m = 0, with m being the mass, and 
another with g^ = (47r)^ and arbitrary m. 

In the following sections we generalize the saddle point equations ( |2.7D to the case 
when the fields p and A are not constant. 



3. Case of non-constant fields - general approach 

In this section we introduce a set of equations that govern evolutionary processes in the 
theory with Euclidean action given in (2.1) in the large limit. These equations are 



2U'{p{x)) =iX{x), p{x) = Greg{x,x), (- + iA(x)) y) = 5(x - y) . (3. 1) 
Here G{x,y) is the Green's function and Greg{x,x) is the regularized Green's function 

Greg{x, x) = lim f y) - ——^ . (3.2) 

x-^y \ Att\x — y\y 

In order to derive these equations recall that, as described in section |2|, any constant 
solution of the theory (|2.l|) in the large N limit is a solution of the gap equations (|2.7| ) 



2U'ip)=iX, p = tr (3.3) 

The trace in the second equation is tr = G{x,x), where G{x,y) is a propagator of a 

scalar field in a background iX which solves the last equation in (|3.1|). The regularization 



required for the divergent G{x,x) leads to (3^). The above equations remain valid also 



for non-constant fields. Any solution of these equations will describe some process in the 
theory. Constant solutions described in section ^ are examples of solutions, but there are 
many others. Some examples will be given in what follows. 

The main difficulty in finding solutions stems from the fact that the last equation in 
( |3.1| ) involves computing a Green's function of the operator with an unknown function 
A(x). Since, to the best of our knowledge, there is no closed expression for such a Green's 
function, we shall consider each case separately. 

Yet another approach is suggested in ||9|, |lO|, where the 0{N) gi<j)^ theory in four 
dimensions is explored in a large-A^ limit. In particular, the effective action and the cor- 
responding gap equations are derived within the subspace of slowly varying functions X{x) 
for which the inequality |(9A)^/A^| << 1 holds. In Appendix^ we build upon this idea to 
construct also an effective action in this approximation and analyze different scenarios of 
dynamical evolution within it. 
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4. Time-dependent solution in the conformal case 



In the first case we consider the effective potential 
dimensionless coefficient, the bare potential is thus 

96 76 



U 



6N^ 



?2|) which contains only - the 



(4.1) 



whereas the full quantum effective potential according to (|2.8|) is 



Ueffip) 



9e-9c ^3 



where, as before, p = {(p'^)/N and gc is the critical value of ge 

9c = (47r)2. 



(4.2) 



(4.3) 



As mentioned in section |2|, the theory is scale invariant also quantum mechanically and 
possesses two different phases: one with unbroken scale invariance for g^ < gc, and the 
other with spontaneously broken scale invariance and a spontaneously generated mass for 
96 = 9c- The behavior of the system depends on the sign of g^ in the classical case and in 
addition on whether g^ is bigger or smaller than gc in the quantum case. Note that from 
( |2.7| ) it follows that relevant values of p after a regularization are negative. The quantum 
potential for both signs of gg — gc and some types of motion in it that will be considered 
below are shown in figure |l|. This potential plays a role in various attempts to obtain 
a holographic dual to a gravitational system which may be suspect of exhibiting a Big 
Crunch 0] . 

The unstable potential resulting for g^ > gc can be misleading. In that case the UV 
cutoff A can not be removed, a UV completion is needed to have a cutoff independent 
theory. In this paper we check what would have been the time scales involved if one 
insisted to keep the unstable potential as is, with neither a cutoff nor a UV completion. 
This is done to learn about the properties of such unstable potentials when they actually 



arise. 





a) 



b) 



Figure 1: a) Effective potential for — gc — I, b) Effective potential for gg 



In this section we present a solution to equations ( |3.lD in the scale invariant phase. 
We consider a situation when all fields depend only on one coordinate out of three, and we 
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choose it to be a time t. The scale invariance of the potential dictates the asymptotic form 
of the time-dependent solution. Initial conditions such as the initial time Iq and the initial 
value of the field at time to break scale invariance and are thus reflected in the solution as 
well. The solution will require two integration constants. 

The solution that we find is actually the classical motion described by the full quantum 
effective action T{(j)). We compare this solution with its classical counterpart, by which 
we mean a classical motion in the quantum effective potential (^]^) with all corrections 
that include derivatives of the field neglected. Since p ^ (j)'^ this last case looks as a purely 
classical problem of the motion in the potential cp^. In what follows we will call this motion 
a classical motion, in contrast with the motion described by the full action T((p), which we 
will call the quantum one. 

Here we study a solution for which the classical energy vanishes. This initial condition 
does not break scale invariance. However, a second initial condition, for example, the value 
of the field cpo at the time to, does introduce a scale. We compute the time evolution of 
the field for both the classical and quantum cases for zero energy and show that although 
qualitatively the two cases are similar there is a quantitative difference: the quantum rolling 
is faster. If one considers a time which takes the field to get from a given initial value to 
infinity in a way that the total energy vanishes then the precise relation between these 
times in the classical and quantum cases is 

A.,^^ ,4.4) 

A solution for the non vanishing energy in the classical case is presented in Appendix |^. 
4.1 The classical case 



We consider here a classical motion of the field in the effective potential ( |4.2| ). We start 
with the following Lagrangian 



L = ^{d,$y-^,m\ (4.5) 



1 

2 V"''^; 6Af2 

where a is an arbitrary constant. For time-dependent and space-independent solutions, the 
Lagrangian is 

where dot means a derivative w.r.t. t. The point cp'^ = in the space of fields corresponds 
to the extremum of the potential. In what follows we mostly consider an inverted conformal 
potential and investigate those classical solutions with vanishing total energy, which either 
asymptotically tend to </>2 = or increase indefinitely. Such solutions are closely related to 
the exact quantum ones which we are going to find in the next subsection, in a sense that 
for both cases (classical and quantum) the total energy is zero, the system approaches the 
equilibrium point c/^^ = as t — > ±00 and the field's magnitude diverges after a finite time. 
The general ansatz for any classical solution of the considered problem is given by 



(l) = VN(t){x)h, (4.7) 
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where </'(x) is an arbitrary function and n is a unit vector in the space of fields. For the class 
of solutions which we are interested in, </>(x) can be determined by the scaling arguments 
(see below). Moreover, one can show that n = const within the subspace of solutions under 
consideration. Substituting the above ansatz into the Lagrangian and bearing in mind that 
n = const, one obtains 



6 



The corresponding EOM is 



(4.9) 



By scaling arguments the solution of vanishing energy should be of the form 



^ pa = ±-^^ , (4.10) 



^/±{t-tdiv) t - tdi 

where t^iv is a time when the field may diverge and pd is an analog of p, which is (f)^ /N . 
For a positive value of the coefficient a of the (j)^ potential the only zero energy solution is 
that in which the field is stuck statically at the minimum of the potential. To obtain truly 
time-dependent solutions with zero energy one needs to study the case of the potential 
unbounded from below obtained for a < 0. In this case the EOM determines the value of 
7 to be 

Also in this case, there are two different kinds of time-dependent solutions with zero energy: 
the field can either escape towards infinity or asymptotically tend to the maximum of the 
potential. We will consider the first case, which implies that in the solution ( [4.10| ) there 
should he t < tdiv and the sign should be minus. Insisting on placing the field at the 
maximum of the potential at time to will result with a static solution remaining there. In 
order to obtain a time-dependent solution we can set either (po or 0o to be non-zero. The 
one determines the other through 

With these initial conditions the solution is 

7 7^ 



Pel = '—-T , (4.13) 



Po ^0 



with 7 as in ( 4.11 ). From here one reads off that 



7^ 

tdiv =to + jo- (4.14) 

That means that the field starting at a finite value (po at time to reaches an infinite value 
after a finite time At^i which is 



1 
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In particular, for the potential (4^) this time is 



At^ = ^J^— . (4.16) 

This case can be interpreted as a calculating the divergence time in an exact quantum 
effective action in which, however, no account was taken yet of the time derivatives of the 
fields beyond the standard kinetic term. We next turn to consider the effect of this time 
dependence. 

4.2 The full quantum solution 

Here we consider the full quantum evolution with all time derivative corrections included. 
As in the previous section, we consider a process with vanishing energy. Similarly to the 
classical case, the scale invariance fixes the time dependence of the functions 

Pit) = , ^X{t) = .2 ■ (4-17) 

where the subscript L stands for "Lorentzian." In this section we compute the coefficients 
(XL and ul- We show that in addition to the trivial solution = cti = there exists in 
the case > Qc another one 

2 Vffe - Qc 4: 96- 9c 

If qq > Qc then the quantum potential is unbounded from below in the sense described 
before, and the nontrivial solution represents a rolling from the maximum of the potential 
towards infinity. However, if < (/g < 5c then the potential has a minimum at p = with 
vanishing energy. Since we consider only solutions with zero energy, we find only a trivial 
one, corresponding to a field at the bottom of the potential well. 

From this one can compute, as for the classical rolling time At^, the time it 

takes the field to roll from a given value of the field /?o to an infinite value, according to 
the solution ( [4.17 ) with coefficients ( 4.1^ ); this time turns out to be 

1 

2po\/96 - 9c 



At, = ^ . (4.19) 



One can compare this rolling time with a similar time Atd for a rolling in the effective 
potential (ff^) with p replaced by cp'^, as explained in the previous section. The result is 



At, 



At, = ^ , (4.20) 

the quantum rolling is faster than the corresponding classical one. It should be stressed 
that there is no reason for the two times to be equal, since the effective potential captures 
only those terms in the full quantum effective action which do not involve time derivatives 
of the field beyond the canonical kinetic term; those terms do however appear in the case 
of non-constant solutions. We will see another example of this phenomenon in the next 
section, where we consider a rolling in a general renormalizable potential. 
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We start all computations in the Euclidean signature. The scale invariant ansatz for 
the functions p and iA is 

pi^r) = , i\{T) = , (4.21) 

where r is a Euclidean time and r^^ji, is a Euclidean time instance at which the field diverges. 
The subscript E emphasizes that the work in the Euclidean space-time. In what follows we 
will choose Tdiv to be zero. This choice fixes a time translation symmetry of the problem. 

The Lorentzian motion is recovered by replacing r by it. This, in turn, leads to the 
expression ([4.17] ) for p with = —iuE- The field p is real if a^; is imaginary. Since, 
as explained in section ^, p is a value of the regularized Green's function with coincident 
points, the full Green's function must be complex. 

The next step is to check if there is a non real-valued Green's function of the operator 
— □ -|- a£;/r^. There are two different cases: 1) ue < —1/4, 2) ue > —1/4. The Green's 
functions for both cases are derived in Appendix]^. 

In the case aE > —1/4 the Green's function is 



where 



= ^^^. (4.23) 

and the Euclidean times r and tq must be of the same sign. This Green's function is real 
and therefore does not lead to a real-valued field p. This case corresponds to the potential 
(p^ with a positive coefficient, which is bounded from below. In this case the only solution 
with zero energy is the trivial one p = 0. 

Turn now to the case < —1/4, which, as will be shown in course of the derivation, 
corresponds to the potential unbounded from below. In this case /3 in ( [4.23 ) becomes 



imaginary and the Green's function ( 4.22| ) becomes complex- valued.^ The short distance 



expansion of this function (for r ^ and r — > tq) is 



G(r, r, to) ^ ^ = - . (4.24) 

47rY^r2 (|r| - |ro|)2 47r|ro| 



The first term is a singularity of precisely the form anticipated in ( |3.2D , and it is canceled 
in the regularized Green's function, which therefore is 



G„,(.,.) = --i^.-,4^. (4.25) 

47r r Svr r 



Since according to (|3.1| ) Greg{'T,T) = p{t) we get 



UE = -I z signr . (4.26) 

ovr 



^As explained in Appendix^, in this case there is also another real-valued Green's function. It describes 
a p = solution in the potential unbounded from below. 



-9- 



The value is imaginary, and therefore after the Wick rotation the solution is real. The sign 
of T which appears in this expression reflects the fact that, as in the classical case, positive 
or negative r's describe the field that either runs to infinity or tends asymptotically to 
zero. From the first equation in (3.1), which in this case is 



ge p\r) = iX{T) , (4.27) 

we conclude that 

aE = \ . (4.28) 

We see that aE > —1/4 holds for ge > gc, and this is precisely the case when the effec- 
tive potential is unbounded from below and one expects to have a rolling solution. The 
coefficient oe defined in ( [4.26 ) is 



aE = - „ , signr . (4.29) 

2V56 - 9c 

Wick rotate now to the Lorentzian signature, t = it. The Lorentzian solution is the 
analytic continuation of our result 



We see that p is negative, as necessary. The evolution of the field is either a rolling from the 
top of the effective potential, which has no minimum in this case, for i < 0, or a running 
towards zero for t > 0. We will concentrate on the case of rolling. Plugging in the value of 
aE from ( [4.28| ) and taking t < we obtain 

p{t) = —^^=, t<0. (4.31) 
2t V56 -9c 

This, after restoring tdiv, leads to the following final form of the rolling solution 

Pit) = , ^m = JT^^ , (4.32) 

where 

(XL = „ , , = 7 ' (4-33) 

2 Vff6 - 9c 4: 96- 9c 

and we have restored the divergence time tdiv 

If 56 is sufficiently close to gc, then all the fields involved in the problem are slowly 
varying functions of time and the approximation suggested in is applicable. We build 
upon this approximation in Appendix ^ in order to validate the exact results obtained 
here. We find a full agreement. 

Moreover, we want to compare these results with the classical ones of the previous 
section. In order to do it we consider the quantum process which is analogous to the 
process considered there: the field starts at certain value po at the initial time to and 
reaches the infinity after some time Atq. This rolling time is 

Atq = — , ^ . (4.34) 

2po V56 - 9c 
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Atq = . (4.35) 



Comparing this rolling time with its classical counterpart ( [4.16| ) we see that the quantum 
rolling is faster 

At,; 

As explained above, the two times need not be equal because in calculating At^ one takes 
into account terms in the effective action that contain derivatives of the field (in addition 
to the usual kinetic term). These terms are not taken into account in the classical rolling. 

In Appendix ^ we show that one can get a feeling whether quantum effects tend 
to accelerate or decelerate a rolling. The leading order in h correction to the force is 
— (T^[/"'(0)/2, where o is the standard deviation. In particular, in our case according to 
equation ( |4.8| ) the correction is — 10acr^(/>^. We see that the classical force —a<\^ and the 
quantum correction are of the same sign, and therefore the quantum correction increases 
the force. So indeed the quantum rolling should be faster. In the computation we have 
determined the precise relation. 

4.3 Energy conservation 

For completeness we illustrate the energy conservation for the process under consideration. 
The energy-momentum tensor in the Euclidean signature is given here by 

r^.(^) = 5;.0a.<?- V[^(5^<?)' + ^(<?')']. (4.36) 

The expectation value of the energy-momentum tensor is calculated by noting that to 
leading order in the expansion 

and that 



< d^^{x)d,^{y) >= Ar^^_^G(x, y) . (4.38) 



Using ( 4.22| ) we get the following regularized expressions 



hm |iG(r,r,ro)= lim l|-G(r, r, tq) = ^^j^-^^, (4.39) 
where 6t = t — tq. 

In what follows we need to compute only the Tqq{x) since we are interested to demon- 
strate only the energy conservation (recall that we consider field configurations which 
depend only on time). First note that due to rotational and translational invariance in the 
r-plane one has 

< mx)d^^iy) > = iV_|L_G(r,r,ro) = -iV^^G(r, r, tq) = -iVAG(r, r, tq) 
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As a result, combining ( [4.21 ), ( 4.3(: )-( 4.41 ) altogether, we finally obtain 



(Too) 



— lim 



(- 



- G{r, r,ro) + {-^ + A) G(r, r, ro)) " f = . (4.42) 



2 r,5T^o ydrdTo 

The zero value of the energy of the system is indeed conserved during the process considered. 



5. Case of small changes of fields 

We next consider a case of renormalizable po- 
tential with non- vanishing coefficients g2, Qa and 
(the restriction to this particular renormalizable po- 
tential is not necessary and can be omitted without 
any loss of generality). For our needs we choose 
52 < 0, 54 < 0, < < Qci where Qc is defined in 
( |4.3| ). The corresponding effective potential accord- 
ing to (|]D is 




92 , 54 2 , f 6 - 5c 3 



(5.1) 



Figure 2: Effective potential (5.1) for 
£,2 = -10, 34 = -100 and ge = 10. The 
field p is negative after the renormaliza- 
tion 



In this expression the field p is constrained to be 
negative. The potential is drawn in figure |^. It has 
a metastable ground state at the origin and stable minimum away from the origin. 

In this section we describe another approach to obtaining time-dependent solutions. 
It is valid in an approximation in which the field does not deviate too much from the 



extremum of the potential. The potential (5.1) possesses two extrema 




1 ± 




(5.2) 



These extrema are shown in figure |2|. 

We study the solution rolling from the top of the potential down to its minimum 
and the solution which oscillates around the minimum. We derive a characteristic time 
of runaway from the top of the potential /?q as well as the frequency of small oscillations 
around the true minimum . These will turn out to be determined by the equations 

LO 47r 



arctan qr = + ■ (5-3) 

SttPq 54 + 256 Po 

In the case of oscillations this equation determines the frequency uj, whereas for the rolling 
it gives the inverse rolling time t ~ 1/uj. Relying on (5.2) one can show explicitly that this 



equation possesses real solutions along with imaginary ones, depending on the sign before 
the square root in ( p.2\ ). We explain their physical meaning below. 

We also compare the results with the corresponding classical counterparts for the same 
effective potential in the approximation when the two extrema are close to each other 



Po - Po 
Pt 



< 1 



(5.4) 
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We find that within this approximation the frequency of the quantum oscillations iVq is 
larger than that of the classical ones tOd 



g 



VSuJci ■ (5.5) 



The characteristic quantum rolling time for tq is smaller than its classical counterpart t^i 

,5.6, 

We turn now to the derivation of these results. For the choice of potential equations 
(|3.1|) are given by 



g2+gip{x)+gQp (x) = i\{x), p{x) = Greg{x,x), y-ax+iX{x) j G{x , y) = 6{x-y) . 

(5.7) 

Let us assume that the total change of p and A in the course of dynamical evolution is 
small compared to their mean values. This is the case if one would explore the system 
during short time intervals. It also holds during the whole evolution process if one imposes 
condition ( |5.4| ). 

We expand the fields as 

5p{x) = p{x) - po = Cpi{x) + C'^p2{x) + 

5\{x) = \{x) -Xo = C Xi{x) + X2{x) + (5.8) 



where po and iXo provide a constant solution to (^J) and C is an expansion parameter 
which is assumed to be small: |C| <C 1. As mentioned in section ||, in the case of a 
constant solution at the minimum of the potential, the value of iX coincides with a square 
of a physical mass, and therefore 

zAo = . (5.9) 



The last equation in (5^) can be written as 

(^-□^ + m^ + i5X{x)'^ G{x, y) = 5{x - y) . (5.10) 
The solution to this equation can be found by performing a series expansion 

G{x,y) = Gmix,y) - j (fzGra{x,z)i5X{z) Gm{z,y)+ 

+ d^zcfwGm{x,z)i6X{z)Gmiz,w)i6X{w)Gm{w,y) + ..., (5.11) 



where Gmix,y) is a propagator of a free massive scalar field 

-m.\x-y\ 

Gm{x,y) = — (5.12) 

4:Tr\x — y\ 

One needs to regularize Green's function at coincident points. In the series expansion 
(5.11) only the first term diverges in the limit x ^ y, and its regularized value is 

/ ^-m\x-y\ ^ \ ^ 
Gm,reg{x,x) = Hui — r - — r = -— . (5.13) 

x^y\ 4:7r\x — y\ 47r|x — y| / 47r 
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This result, in particular, shows that po = —m/4:ir in accordance with the known result 
presented in section |2|. Altogether, the expression for p becomes 

p{x) = —- / (fizG'^{x,z)i6X{z)+ 



Air 



+ / (fz(fwGm{x,z)i6X{z)Gm{z,w)i6X{w)Gmiw,x) + ... (5.14) 



One can solve equations (5.7) order by order in C. 
• Order 0. To this order equations ( |5.7D yield 

2 2 

92 + 94Po + 96Po = m , = - — , (5.15) 

where, as above, zAq = m?- There are two solutions: 



^° = 2(/-ff6) (^ ^ r + ^ " 1 ' ^^-^^^ 



Order 1. In this order one has to compute the integral in the second term of the RHS 
in ( ^.14 ). It involves a convolution of and zAi, and the result is 

1 

Piix) = ; arctan iXAx) . (5.17) 



The first equation in (5^) then gives 

\ 47r^/— □ 

54 + 2po5'6 ) Pi (a;) = 7^ Pi (a;) . (5.18) 

^ arctan ^ — 

If we look for a solution which depends on a single coordinate, a Euclidean time r, 
this equation simplifies 

)47r zc^ 
arctan^ 

and its form suggests the following solution 

Pi(r) = e*-^ , (5.20) 

where w satisfies equation ( |5.3D . 

The result for oo depends on the choice of po^ ^ is real for pQ = p^ and imaginary for 
Po = Pq ■ This means the following: p^ corresponds to a maximum of the potential 
and the result describes a runaway solution, whereas p^ is a minimum of the potential 
and in that case u; describes the frequency of the oscillations around it (recall that 
the computation was carried out in a Euclidean time and in order to get a Lorentzian 
solution one has to Wick rotate the time, which introduces an additional factor i in 
the exponents). 

The solution for iX to this order is 

iXi{T) = Aie'^\ Ai= ^^"^ ^ =g^ + 2ff6Po- (5.21) 
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One can continue the perturbative expansion to higher orders. There will appear 
corrections both to the frequency of oscillations and to the shape of the solution. This 
follows from the fact that if the amplitude is large enough the field will overshoot the value 
Pq and there will be no oscillations around it. 

Let us look more closely at the frequency to of small oscillations around Pq . It is 
determined by equation ( |5.3D . Assume that the two extrema of the potential are close to 
each other, as in equation (p-4i). This approximation will be valid if the following condition 
holds 

^(5c-<76)>-l. (5.22) 

at 

Due to the assumption about the coupling constants (52 < 0, 54 < 0, < (/e < dc), this 
criterion can be satisfied. In this case to is small and thus it can be determined by a simpler 
equation obtained by expanding the arctangent to the third order 




(5.23) 

This is the Euclidean result, the Lorentzian one will be of opposite sign. The Euclidean 
result for ui^ is negative for qq < g^, and therefore the Lorentzian result will be positive, as 
expected. In Appendix]^ we rederive this result in the framework of another approximation 
presented in [^. We find full agreement between the results. 

This frequency can be compared to a frequency iVd of oscillations of a classical field 



(p in an effective potential (5J.). Since the potential is 0{N) invariant the most rapid 
oscillations occur in the radial direction, thus the maximal value of the classical frequency 

is determined by the equation o;^; = U"{(I)q), where (po = \I~P^ is the value of the field (p at 
the extremum. In our approximation the Euclidean result is 

^l^ ^./l + %(5c-56). (5.24) 

Therefore 

~ V^UJcl , (5.25) 

and we conclude that the quantum field oscillates faster around the true minimum than 
the classical one. 

One can carry out a similar computation concerning the characteristic time of escape 
from the potential maximum at p^ . The result is similar: the escape time in the quantum 
case is shorter than its classical counterpart by the same factor of \/3- 

6. Tunneling 

In this section we investigate a possible tunneling in the system. Following Coleman we 
search for solutions which depend on a Euclidean radius 



rE 



2 + r2, (6.1) 
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where r is a distance in space and r is a Euclidean time. The physical meaning of these 
solutions is revealed after the Wick rotation back to the Lorentzian time. The rotated 
solution depends on a combination = 2 — where t is a Lorentzian time. Such 
a solution describes a spherical wave in the space time. Within the semiclassical picture 
of the field evolution after the tunneling (an expanding bubble with a thin wall) the real 
and positive values of correspond to points where the wall has not yet arrived, whereas 
imaginary values correspond to those points where it has already passed. 

In what follows we will work in the Euclidean signature and will suppress the subscript 

E. 

Relying on the aforementioned feature |^] that three dimensional theory in the large 
N limit is renormalizable regardless the shape of the potential, we base our approach on 
the reverse engineering method: we will derive the equation for the instanton, then we will 
choose a solution which describes a tunneling in some unknown effective potential, and at 
the end we will reconstruct numerically the potential itself. Our consideration is reliable if 
the false and the true vacua are close to each other. 

Our main results are summarized in figure ^, where figure |^ is a plot of a function 
which describes a deviation of the solution from its value at the false vacuum as a function 
of the Lorentzian radius r = yjr'^ — t'^ for its real values (r = 00 corresponds to the false 
vacuum itself) , and figure §3 is a plot of the effective potential with false and true vacua. 
We also compute the tunneling amplitude and compare it to the corresponding semiclassical 



<5p Ueff 




a) b) 



Figure 3: a) Function bp[r) giving a deviation from the false vacuum which is at r = 00, b) 
Resulting effective potential C4//- 

value. In both cases the amplitude is A ~ exp{—NS). In the semiclassical case Ssd — 4.7, 
whereas in the full quantum computation S ~ 3.56, so the quantum amplitude turns out 
to be larger. 

We consider a situation when the change in fields during the tunneling process is small 



Pt-Pf 



Pt 



« 1 , (6.2) 



where pt and pf are the values of p in the true and false vacuum respectively. We start 
with an arbitrary potential U{p) and carry out the solution which presents a combined 
numerical description of a tunneling and a subsequent rolling. The numerical analysis is 
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done relying on the equations which are vahd up to a leading order in the aforementioned 
small parameter. 



Write similarly to ^5. 



~ -— + Cpi(x), iX{x) = m"^ + C iXi{x) . (6.3) 



Air 

To order we will have similarly to ( 5.15| ) 



2U'{po)=m^ Po = -^- (6-4) 



The order 1 equation ( 5.17| ) also remains intact. The two equations can be combined into 
a single one 

2U'{p{x)) = -{ATrfpl + ^^^^ p(x) . (6.5) 

arctan ^ 

This is a general form of the equation which describes a behavior of the system to the 
leading non-trivial order. In the case that we want to consider, namely, when the field p 
depends only on the Euclidean radius r the equation simplifies. As we show in Appendix 
0, on functions depending only on the radius, the operator ^— □ becomes 

^/^p{r) = -Ap{r) = - — PV [ i-^^fli-ds, (6.6) 

vrr J s — r 

— oo 

where PV denotes the Cauchy principal value and we assume that the function p{r) has 
been continued to the region of negative r as an even function 

p{-r) = p{r) . (6.7) 

Equation ( |6.6D also serves as a definition of an operator A which will be extensively used 
in what follows. 

From what was said above we arrive at the final form of the equation for p{r) 

Att a 

2U'{p{r)) = -{A^fpl + J- p{r) . (6.8) 

arctan ^ — 

Svrpo 

This equation should be supplemented by boundary conditions for p{r). Following |Q] we 
require that p{r) go its value at the false vacuum when r — > oo. We denote this value of p 
by po and define 

5p{r) = p{r) - Po . (6.9) 



The function 6p{r) goes to as r goes to infinity. With this definition equation (p^ ) 
becomes 

Att A 

2U'{po + 6p{r)) - 2U'{po) = j- Spir), 2U'{po) = (Anfpl . (6.10) 

arctan 
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Linearized around po 

IttA 

U"{po) 5p{r) = J- Sp{r) . (6.11) 

arctan „ 

Note, that since at po the potential has a minimum U"{pq) is positive. However, as shown in 
Appendix]^, all eigenvalues of A are negative, and, since po < 0, it follows that the operator 
on RHS of (6.11) is negatively definite. Therefore there are no solutions to equation (6.11). 
Nevertheless, this feature may not pose a problem since one cannot extrapolate it to the 
solution of equation ( |6.10| ). 

One can rewrite equation (6.10) in the following form 

A-TT A 

U"{po)6p{r)+a{6p{r)) =B5p{r), B = ^, (6.12) 

arctan ^ — 

where the function cr((5/9(r)) is proportional to 5p^{r). Since for large real r the function 
5p{r) goes to we conclude that it should become an eigenfunction of B at least for very 
large r. 

We say that a function /(r) is an "asymptotic eigenfunction" of the operator B if there 
is a number k (an "asymptotic eigenvalue") such that / obeys the following requirement 

hm /(r) = 0, hm ^ f i"^) - ^ f = q . (6.13) 

r— >oo r— >oo f{f) 

In Appendix ^ we show that the operator 13 possesses many asymptotic eigenfunctions 
with positive asymptotic eigenvalues (see ( [E.15 ) and a discussion there and remember that 
Po < 0). Consider the following illustrative example. 

The construction of the candidate function 6p starts with a function /i(r) 

L 0, r > 1 

This function is plotted in figure Without any loss of generality, let us take both the 
eigenvalue and po to be equal —1. Then one can compute (see Appendix^ for details) the 
functions 6p and I36p 

oo 

^ , , 2 f ,,, arctan^ sinkr „ 

5p{r) = / ^{k)— dk , 

TT J Airk + arctan ^ r 


oo 

A,- / N 2 /" , Ank sinfcr 

B6p{r) = - / ^{k)— jr dk , 6.15 

TT J Airk + arctan ^ r 



where ^{k) is the sine-Fourier transform of /i(r) 

, , , , , 40A;(/c2 - 15) cos /c - 120(2A:2 - 5) sin /c 
^{k) = (r) sm{kr)dr = ^ ^ ^ ^ . (6.16) 
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These functions are plotted in figure and |b. One can also reconstruct the classical 
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Figure 4: a) Function fi{r), as defined in (6.14), b) Function B5p{r). 



potential indeed increases from its value at p = 1 (the fact that its derivative does not 
vanish at this point is a numerical error, the numerics is not reliable at this region; the 
derivative is very small though). We see also that the effective potential has a minimum 
at yO ~ 1.012 and then grows up. Since p — po <ti p the approximation is valid along the 
process . 

These results have the following interpretation. At t = the field tunnels from the 
false vacuum and a bubble is created. Its shape is given by the function p{r) that we 
computed, and there are points in space where the field acquires all values that are covered 
by the plot of the potential Then, as the time grows, the solution evolves according to 
the law p{V — t'^). We see that the initial shape of the bubble is complicated, contrary 
to the case of the thin wall approximation. 

The fact that our numerical computation reconstructs the minimum of the potential 
and even goes further means the following. In the semiclassical approximation the compu- 
tation similar to the one presented here would reconstruct the potential up to a point which 
is just a little beyond the point where the value of the potential is equal to its value at the 
false vacuum. That the field at all gets beyond the turning point (a point where the poten- 
tial is precisely equal to its value at the false vacuum) is because of the friction term in the 
EOM. If there is no friction term, like in the case of quantum mechanics, the computation 
will only be able to reconstruct the potential up to its turning point. Our computation 
gets so far beyond the turning point and even beyond the true vacuum because the action 
we work with is very complicated. In particular, it possesses a non-standard kinetic term, 
which, in turn, leads to a non-trivial friction.'^ 

^We are grateful to J. Barbon for discussions on these issues. 
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Having computed the potential, one can also evaluate the effective action on the solu- 
tion in order to get the tunneling amplitude. In the case under consideration the action is 
Seff ~ 3.56, and the tunneling amplitude is exp{—NSeff)- 

This result can be compared with a prediction of the semiclassical approximation for 
the same potential. In that computation one has to solve a radial Euclidean EOM for a 
bounce and then to evaluate its action. The numerical result in this case is Ssc — 4.7, and 
the tunneling amplitude, similarly to the case above, is exp(— S'g//). We see that the 
full quantum tunneling amplitude is larger in this case. Once again this falls within the 
pattern that taking account of time-dependent effects shortens the characteristic time of 
the process. 

7. Conclusions 

We have studied a variety of time-dependent processes. These included rolling of fields, 
tunneling among vacua and their decay. This was done for a large class of vector 0{N) 
models in three space time dimensions, among them the system which is conformal also 
quantum mechanically. Using the methods of the large N expansion we were able to take 
into account the effects of the time variation of the fields on the exact effective potentials. 
We have calculated exact and approximate characteristic time scales of such processes. In 
most cases the results were quantitatively different than those obtained without considering 
the time dependency in the effective potential. A qualitative difference is found when 
analyzing the bubble driving the decay of a false vacuum. The bubble shows thick rather 
than thin characteristics. For the cases studied a pattern of accelerated time scales emerged. 
The next natural step is to study supersymmetric extensions of such systems and their 
coupling to gravity. 
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Appendices 

A. Effective action in the region of slowly varying fields 

In this appendix we follow |^] in order to compute the effective action which governs the 
dynamics of the large- theory within the subspace of slowly varying functions A(x), for 
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which the inequahty \ {dX)'^ / \'^\ << 1 holds. We compare the results obtained for the exact 
conformal case and those obtained for small field variations. The comparison is done in 
those regions where the validity of the approximations overlap. 

For simplicity of notation we suppress imaginary unit i in front of A in what follows, 
that is iX — A. 

Let us define r(A) from the Euclidean functional integral 



=-r(A) 



N 



1 



-□ + A)(/) 



exp 



Trln(-n + A) 



(A.l) 



In the regime of aforementioned approximation r(A) can be written as a local expansion 

r(A) = j d\ [Fo(A) + F^{X){d^Xf + . . .] , (A.2) 

where -Fo(A) and Fi{X) are local functions of X{x). In particular, Fo(A) is found by calcu- 
lating the Tr ln(—n + A) with a constant A. Up to an infinite constant, which is zero in 
the dimensional regularization, we have 



and thus 

W) = -— A3/2. (A.4) 
On the other hand, in order to calculate -Pi (A) one notes that 



,5^r(A) 

'a(x)5A(y) 



1 

4 
iV 

y 

N 
IGtt 



(27r)3 

d^k 



d^p 



[All 
J W 



(27r)3((p-A:)2 + A)(p2 + A) 

''^"^-y^ [' da {a{l - a)k^ + Xy'^^ 
Jo 



(A.5) 



This can be expanded in the form 



S^TjX) 
'SX{x)6X{y) 



+ 0{n^6{x - y)) , 



(A.6) 
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where in the second equahty we assume that A is constant. Such an expansion is justified 
due to the assumed approximation regime of slowly varying fields. Using equation (A. 2) 
we can now identify Fq{X) and -Fi(A) with the constant-A expansion of equation ( [A.£| ) 



92 A 
Fi(A) 



N 
N 



^-1/2 



3847r 



A 



-^3/2 



(A.7) 



Integrating the first of these equations reproduces the result of ( [A.^ ). Thus, effective action 
( |2.6| ) can be rewritten as follows 



Seff{p,X) = d X 



2 127r 3847r ^ ^ 



If we now define 



|l/4 



then effective action can be written as follows 

(12^)2 



Seffip,'4') 



d^x 



u{p) - '-^Pi^' - (i2^)V' + i^m'' 



(A.9) 



(A.IO) 



Let us exploit this action in order to illustrate the process of tunneling and rolling in 
the case when the potential is given by 



52 , ff4 2 



(A.ll) 



Under this assumption the effective action S(.ff{p^ ip) turns out to be quadratic in auxiliary 
field p and thus integrating it out yields 



S^ffi:^) = [ d'x \\{di.f - ^^V'S - (12vr)V^ + |^(12vr)V^ 

J [2 4^4 2g'4 



(A.12) 



Let us explore the temporal rolling first, that is we consider the situation when all the 
fields are time-dependent only. As a result, the corresponding Lorentzian equation of 
motion represents a particle of unit mass moving in a potential 



4c/4 2^4 



Thus, 



E 



1 



2\ dt 



(A.13) 



(A.14) 



is a constant of the motion. This can be used to determine the qualitative features of the 
solutions by inspection. 
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As a simple example of rolling, consider the potential shown in figure || which corre- 
sponds to 

52 < 0, 54 <0. (A.15) 

We are interested to investigate the rolling of the system from the local maximum of the 
potential situated at 

9a 



967r 



SvrV 



down to the false vacuum located at 

54 



967r2 



(A.16) 



(A.17) 



The corresponding runaway and oscillation fre- 
quencies are given respectively by 



v[iA] 



(A.18) 



In this case, we are dealing with a solution of the 
equation of motion with E = V{^l^-), whence ^/J as 
a function of t is given implicitly by 
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Figure 5: The effective potential ( A. 13 ) 
(A.19) for 54 = -(12^)^ 52 = 0.91(54/8^)^. 



where t+ is an integration constant, the time at which ■0 equals 

On the other hand, according to Coleman if we are interested to compute a decay 
probability per unit time per unit volume T /V , of the unstable state '0+ due to the barrier 
penetration, one must find the bounce '0, a solution of the Euclidean equations of motion 



ar^ r ar 

subject to the following boundary conditions 

lim ^/5(r) = , '0'(r = 0) 



0. 



(A.20) 



(A.21) 



where prime denotes derivative with respect to Euclidean radius r. To leading order in 



T/V 



-N 



s-ffi^){l + 0{l/N)) 



(A.22) 



In other words, in the limit of large N barrier penetration is exponentially small and 
thus, as emphasized in the text, the dynamics of the system is governed by the rolling 
processes only. 

It is instructive to perform the calculation in the region where the approximation of 
slowly varying fields presented above overlaps with the small field approximation used in 
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the section Therefore let us assume that condition ( ^.22| ) holds. As a result, the two 
extrema of the potential plotted on figure ^ can be made arbitrary close to each other, 
and the field derivatives in turn, during the evolution of the system from the top of the 
potential down to its minimum, become arbitrary small. 

As an illustrative example, let us derive the frequency of small oscillations around the 
true minimum . Relying on ( A.10| ) the corresponding EOM are given by 



□-0 



-2(127r)2pV^ - 6(127r) V^ 

„2 



92 + 9iP + gep ■ 

Since the oscillations are small, we linearize around the true vacuum 
6p 



(A.23) 



-2{12TTf{3p^ij^^6tp + i^^^Sp) - 30(127r)2?/;+45V' 



94 + 2gep2 



(A.24) 



where 



3 ' 

Sp = p- p+ ^ 

Sip = i/; — ip^ . 

As a result, the quantum frequency of oscillations is given by 



1 



54 + 256^0 



9c - 96 



1 + 



4^2 

9l 



{9c - 9&) 



(A.25) 



(A.26) 



This result agrees with ( ^.23 ) obtained via small field approximation. 

The other case we will consider is that in which only g^, the dimensionless coefficient, 
is present. Equations ( A.23| ) in this case combine together and yield 



1 




(A.27) 



-6(127r)2 

Since in the region of slowly varying fields 1 >> |((?A)^/A^| ~ |n'i/;/'0^| one concludes that 

89 = 96-9c«l (A.28) 
in order to justify the approximation. This yields 

(A.29) 



uip 



Therefore according to ( 4.10| ) the time-dependent and space- independent solution with 
vanishing energy is given by 



(A.30) 
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with^ a = — 27 6g. 
Altogether 



2„/.4 



A = (12^)^' 



96 



4(56 - 9c) (t - tdivY 
which coincides with the exact results ( 4.171 ), (|4.18D . 



(A.31) 



B. Classical motion in the cj)^ potential 

In this appendix we present a solution which describes the classical evolution of a particle 
with energy E in the 4)^ potential. We derive the divergence time as a function of the 
energy of the particle. 

Consider the following classical Lagrangian 



1-2 o c 
L = -x^ + -x^ , 
2 6 

which describes a run-away potential (for positive a). A general solution is 



^(t) = (9V3-15)^/^ 1^- 



' l-cn(aV6^i/3^,A;) 
2- ^/3 + cn(aV6Ei/3^,A;) ' 



where 



e = 2^/33Vi2(t_to), fc = 1(2 + ^3) 



(B.l) 



(B.2) 



(B.3) 



and to and E are integration constants. E is the total energy and to is an arbitrary time. 




Figure 6: Solutions x{t) for a = 1 and energies E = ±0.25. 

cn(a;, k) is the elliptic cosine function. In figure |6| there are the solutions for a = 1 and 
E = ±0.25. We see that the solution with positive energy passes through 0, whereas the 
one with negative energy never reaches (a particle comes from infinity and bounces back) . 

If a particle starts at some point xq at time to with energy E then it gets to infinity 
after a time 



At 



1 



24/331/12^1/6^1/3 



cn 



cn"^(\/3- 2,/c)- 

1 ^3(3^/3 - 5y/^ E^^ ± 21/331/6(3 _ 2^3) a^/^ : 
3(3^/3 - 5)1/3 ^1/3 + 21/332/3 ^1/3 3.2 



, (B.4) 



*Note that (4.1C) corresponds to the solution of the Lorentzian EOM, whereas (A.23) corresponds to 
the Euclidean time. This reveals the origin of an extra minus sign in the expression for a. 
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where cn~^(x,A;) is the inverse elliptic cosine function. This divergence time is plotted as 
a function of energy in figure |^. We see that the more is the energy the less is the time 
that takes the particle to get to infinity. For small energies this divergence time is 



At 



1 



E 



2x1 V a °'^^"T8fl3/2 



and for vanishing energy one recovers ( 4.15| ) 



(B.5) 



C. Green's function appearing in a conformal 
case 

In this appendix we compute a regularized Green's func- 
tion of the operator 

L = -□ + a/T^ . 

This Green's function is used in section ^. For this we 
derive a short-distance expansion of the Green's function 
which solves the equation 

LxG{x,y) = 6{x - y) . 
The Green's function is given by the following expression: 

^n{x)^*n{y) 




G{x,y)=Y, 



An 



Figure 7: Divergence time as a 
function of energy for a = 1 and 
xo = 1. 

(C.l) 



(C.2) 



where ^''s are eigenfunctions of L and A's are corresponding eigenvectors. 

A Laplacian in L is three-dimensional and involves a Euclidean time r and two more 
coordinates which we will denote collectively by a vector r. L possesses a translational 
invariance in the plane of r and therefore we are free to put the (5-function at any point 
in it, and we will choose this point to be the origin. It also occurs at some time tq. 
With this choice, the Green's function will depend only on the distance from the origin: 
G{x,y) = G(r,r, To), and only angular-independent eigenfunctions will contribute to the 
sum in ( |C.2D . So, the relevant eigenfunctions of L are of the form ^(x) = R{r)T{T). Here 
R{r) satisfies the equation 



-R"(r) 



1 



Wir) = XrR(r) 



Normalizable solutions of this equation exist only for positive and are given by 

R{r)= ^ 
The function T(r) satisfies the equation 



(C.3) 



(C.4) 



a 



T"{r) + -.T(r) = A.T(r) 



(C.5) 
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This equation was investigated in the context of Quantum Mechanics (with r being a 
spacial coordinate) and is known to possess very different kinds of solutions for different 
values of a. We consider first the case of a > —1/4. In this case A,- must be positive and 
the normalizable solutions are 



T{T) = \l'-fJp{^JXr\T\), 



/3 



(C.6) 



In this case there are two disconnected Hilbert spaces for different signs of r (since, in the 
language of Quantum Mechanics, there is no tunneling from a region of r > to a region 
of r < and back), and therefore the full propagator G will vanish unless r and tq are of 
the same sign, and this is what we assume from now on. 
The full normalized eigenfunction of L is therefore 



(C.7) 




^Xr,Xr{r,T) = ^1-1 J^(v^A,|t|) Jo(VArr) , 
and the corresponding eigenvalue is 

A = A^' -\- Xj- . 

According to equation (p.2D, the Green's function is given by 



G'(r,r,ro) 



Using the integrals [|12| 



Svr J A,. + At- 





JoWKr) JpivXr |r|) J/3(VA.r \to 



dx 
x + y 



(C.8) 



(C.9) 



(C.IO) 



dx Ko{a x)Jy{bx)Ju{cx) 



(r+ — r^Y[r^ + r_) 



r± = ^a2 + (6±c)2, (C.ll) 



we end up with the following Green's function of the operator L 



G'(r,r,ro) 



W^-W. 



/3 



W± = Vr2 + (|r|±|ro|)2 . (C.12) 



This Green's function is real if (3 is real, which means, for a > — 1/4. If, however, a < —1/4 
then P is imaginary, and this Green's function is complex. This complex-valued Green's 
function is the one that describes a rolling in the (p^ potential. 

One can build a real Green's function in this case as well. In order to do this we use 
a fact that there is another solution to the Green's function equation ( |C.1| ), namely 



^(r, r. To) 



Vkl |to| 



1 



W+ + W- 



2^ W+ W- \W+- W 



/3 



(C.13) 



-27- 



This function is not appropriate to be a Green's function in the case of real and positive (3 
since it diverges for small Euclidean times. However, for imaginary (3 this function is fine, 
and we can take a linear combination of it and ( |C.12[ ) with equal coefficients to guarantee 
a real result. The sum of coefficients must be 1 for it to indeed be a solution to the Green's 
function equation. So, in the case of imaginary /3 we get the following real-valued Green's 
function 

This Green's function gives a zero solution in an unbounded from below (p^ potential. 
D. Rolling in quantum mechanics 

In this appendix we build on the Ehrenfest theorem in order to compute the leading order 
quantum mechanical correction to the Newton's equations of motion. This in turn helps one 
to conclude whether quantum effects tend to accelerate or decelerate the classical rolling. 

Ehrenfest theorem reveals a quantum mechanical generalization of Newton's second 
law. In particular, for a quantum particle moving in a one-dimensional potential V{x) it 
states 

where m is the mass of the particle and {x{t)) is an expectation value of the position 
operator x{t). 
Let us define 

xit) = , 

m = x{t) - x{t) , 

(x) = —{x{t) +ri) = — + —^rj + -—^if + . . . (D.2) 

ax ax dx ax^ 2 dx^ 



then we get 



Taking the expectation value of both sides yields 

d'^ , . , d ,,,,,, dV Id^V 2 s 

HiL \AiJU lAitJb ^ xJjiAy 

where a = \/Tfp^ is the standard deviation. 

Thus, in order to figure out whether quantum mechanical corrections slow down or 
accelerate the rolling, one has to fix the relative sign between the expressions on the right 
hand side of the last identity. It turns out that this sign is not invariant and changes with 
the shape of potential. 



E. Operator -\/— □ on radial functions 

In this appendix we derive a form of the operator -y/— □ when it acts on functions which 
depend only on a radius and discuss some related questions relevant to section |6|. 
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If a function p{x) depends only on the radius, then the operator simpUfies. 
Indeed, consider the Fourier transform of p{x). It depends only on the absolute value of 
the momentum p and therefore p{x) in this case can be written as 



oo 



f (rp _ . 1 /" 



where in the last equality we carried out the angular integration. The operator — □ 
multiplies each Fourier mode by \p\ , so it acts on p{x) as 



oo 



-n p(x) = — ^ / dpp^ p{p)smpx . (E.2) 

ZTT^X J 







Recall the definition of the Hilbert transform |11] 

oo 

f{x) ^ H[f]{x) = -PV [ ^dy , (E.3) 
TT J y - x 

— oo 

where PV stands for the Cauchy principal value of the integral. One of its main features 
is that it transforms s'mpx to sign (p) cos and cospx to — sin |p| x. Using the Hilbert 
transform, the operator \/— □ when it acts on a function that depends only on a radius 
can be written as 

V^p{x) = -^PV f^-^dy, (E.4) 
TTx J y — X 

— oo 

where we assumed that the function p{x) is defined for negative x to be 

p{-x) = p{x), (E.5) 



as is suggested also by equation 

Define, as in the body of the paper, an operator A to be 



Apix) = ^PV [^yP^dy. (E.6) 
TTX J y — X 

— oo 

Find its eigenvalues and eigenfunctions. In order to do that recall the following features of 
the Hilbert transform: 1) its square is —1 

Hy]{x)^H[H[f]]{x) = -f{x), (E.7) 

and 2) it commutes with a derivative 

H[f]{x) = H[f]'{x). (E.8) 

Now consider the equation for eigenfunctions of our operator A, which can be written as 

Afkix) = -H[{xfky]{x) = kfkix), (E.9) 

X 
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xfkix)) +k^fk{x) = 0. (E.IO) 



where k is an eigenvalue and fkix) is the corresponding eigenfunction. Acting on both sides 
with the same hnear operator once again and using the properties of the Hilbert transform 
mentioned above, we arrive at the following differential equation 

^< 

The solutions of this equation are e^^^^ /x. Since only functions which are finite at x = ±00 
are considered, k is real. By substituting these eigenfunctions into the equation ( |E.9| ) we 
find that the corresponding eigenvalue is negative and equals to — Thus one concludes 
that A is negative definite. The eigenfunctions that possess definite parity are 

(1) sinkx (2), . coskx 

fk (^) = fk W = ■ 

(2) 

Since we assume that p{r) is even one has to drop off (x). 

Next, we consider the operator B defined similarly to equation (6.12) 

B = — ^ . (E.12) 
arctan A 

This operator is positive definite, as A is negative definite, and can be written as 

B = -Bx, B = -, (E.13) 

x arctanh a 

where d means the derivative w.r.t. x. 

Our next goal is to both show that operator B possesses asymptotic eigenfunctions 
with negative eigenvalues, where the definition of such a function is given in ( 6.13| ), and to 
present a general procedure how to construct examples of them. 

If f{x) is an asymptotic eigenfunction of B with asymptotic eigenvalue k, then g(x) := 
X f{x) obeys the following equation 

hm ^9{x)-kg{x) ^ ^ 

Moreover, the requirement that f{x) tends to zero at infinity can be replaced by the 
equivalent requirement that g{x) defined above diverges at infinity slower than x. 

Let us demonstrate how to construct the function g[x) with the above properties and 
negative k. For simplicity of notation we take k = —1. By the definition of g{x) one 
concludes that the function defined by 

h{x) = {B + l)g{x), k<0 (E.15) 



might tend to zero at infinity faster than g{x). From equation ( [E.5D it follows that both 
functions g{x) and have to be odd. Thus, taking an odd function fi{x) with compact 
support guarantees ( E.14[ ) and the following relation holds^ 

00 

g{x) = {B + l)-^h{x) = lj nk)^^^^^^^s^krdk, (E.16) 



^Note that the inverse operator (_B + 1) ^ exists since as we have shown B is positive definite. 
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where ^(k) is the sine-Fourier transform of fi{r) 



fi (r) sm(kr)dr 



(E.17) 
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